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Abstract 



Constrained quantum dynamics is used to propose a nonlinear dy- 



namical equation for pure states of a generalized coarse-grained sys- 



tem. The relevant constraint is given either by the generalized purity 
or by the generalized invariant fluctuation, and the coarse-grained 
pure states correspond to the generalized coherent i.e. generalized 
nonentangled states. Open system model of the coarse-graining is 
discussed. It is shown that in this model and in the weak coupling 
limit the constrained dynamical equations coincide with an equation 
{vq ■ for pointer states, based on Hilbert-Schmidt distance, that was previ- 

iX . ously suggested in the context of the decoherence theory. 

o 
o 

1 Introduction 

^ ' Coarse-grained description of a dynamical system is based on a separation 

of observables into two classes: the class of distinguished i.e. important 
observables and the class of observables that are considered inaccessible. 
A new system, whose state can be maximally determined by the preferred 
observables, is then defined. Evolution of the new, coarse-grained, system 
should be completely described in terms of the distinguished observables only. 
From an operational point of view, the choice of distinguished observables 
is dictated in practice by what can be measured on the given system. For 
example: a) If all observables of a quantum system with the Hilbert space 
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with N complex dimensions H N are considered experimentally accessible 
then every Hermitian operator represents an observable of the system, i.e. 
all observables are distinguished, and the pure states of the system are by def- 
inition rays in the Hilbert space H ; b) In the case of two spatially separated 
qubits H 4 = H 2 <g> H 2 one could consider experimentally accessible only the 
local observables a] <g) 1 and 1 <g> a 2 . In this case the coarse-grained states are 
the product states; c) In a collection of n spins H 2 " = H 2 <g> H 2 . . . H 2 one 
might be able to observe only the macroscopic magnetizations rrii = J2 (T%jn. 
which are then the coarse-grained distinguished observables 

The set of distinguished observables as a subset of the algebra U(N) is 
important in the definitions of notions such as quantum degrees of freedom 
and quantum integrability [TJ, generalized coherent states[21 E] and gener- 
alized entanglement [U El El E|> and provides a framework to study the re- 
lations between these notions [8]. In particular we shall be interested in 
the coarse-grained states representing the generalized non-entangled states 
as introduced and studied in [U El E]- Our goal is to derive an evolution 
equation for which the set of g-nonentangled states is invariant for arbitrary 
Hamiltonian, and discuss its physical interpretation. 

In the next section we recapitulate the theory of generalized entanglement 
and generalized purity. In section 3 we treat a quantum dynamical system 
on H N as a classical Hamiltonian system on R 2N , which enables us to discuss 
constrained quantum dynamics. This is used to derive an evolution equation 
of the states which preserves the maximal generalized purity, that is of the 
g-nonentangled states. This evolution equation is nonlinear and can gener- 
ate, depending on the Hamiltonian, chaotic dynamics of the coarse-grained 
system. Open quantum system model of the evolution of the distinguished 
states is discussed in section 4. In this section we show that our constrained 
evolution equation coincides in the weak coupling limit with the approximate 
evolution equation of the robust states derived in [9] by different means and 
in the context of decoherence theory. 

2 Generalized entanglement and generalized 
purity 

A selected set of distinguished observables is used to define the generalized 
notions of non-entangled and entangled states. The coarse-graining by the 
distinguished observables, understood in the traditional probabilistic sens as 
replacing probabilities by conditional probabilities is crucial in this definition. 
Consider a subset g e u(N) of distinguished observables. A state p g is 



called g-reduced state of the state p if Tr[pL[] = Tr[p g Li] for any L\ G g. The 
reduced state p g is the projection of the state p on the subspace determined by 
distinguished observables. Identifying the quantum states with probabilities 
the standard definition of the conditional probability is recognized. Pure 
state p = \if) >< ip\ is generalized non-entangled if the corresponding reduced 
state p g is pure p 2 g — p g . Otherwise the pure state \ip > is g-entangled. In the 
case that the Hubert space has the bipartite tensor product structure and 
each distinguished observable act nontrivially only in one of the components, 
the previous definition of g-entanglement reduces to the standard definition 
of the bipartite entanglement for pure states. 

In a large class of situations of physical interest the set of distinguished 
observables forms a Lie algebra. In this case a measure of the generalized 
entanglement of the pure state \ip > is provided by the generalized purity, 
which is the purity of the reduced state p g , and is given by: 

^W=E<^I^> 2 , In eg (1) 

i 

where L\ form a bases of the Lie algebra g. The state \ip > is generalized 
non-entangled if P g (ip) is maximal. Pure states with P g (ip) less then maxi- 
mal represent g-entangled states, i.e. the states in which the coarse grained 
system g is entangled with the environment, i.e. with the operators not in g. 
Obviously, whether a pure state \ip >G H N is generalized entangled or not de- 
pends on the choice of the distinguished observables. Once the distinguished 
observables are chosen, the question if the future orbit of a g-nonentangled 
\%f) > will remain in the set of g-nonentangled states depends on the evolution 
equation satisfied by \ip >. 

An equivalent measures of g-entanglement is given by the total dispersion 
of the algebra of distinguished observables 

a,(vo = E < if > - < u > 2 = E(^) 2 - (2) 

i i 

A s (/0) is minimal iff P g {i/j) is maximal. Expressions for the minimal value of 
A s (/0) and the maximal value of P g (ip) in terms of the simple roots of g are 
known [101 [7] and read 

A- g (ip) > E ^ < a '' ai > — ™ n ! Pqi^) — < C 2 > —min = max (3) 

i 

where the highest weight vector X = J2i ki a i i n terms of simple roots o>i and 
C is the quadratic Casimir operator. 

Generalized coherent states have been defined for an arbitrary semi-simple 
Lie algebra. If the algebra of distinguished observables g is semi-simple the 



minimum of A g (/0) and the maximum of P q (ijj) is achieved on the corre- 
sponding generalized coherent states [ 1 1 J . Thus, in this case the class of 
g-nonentangled and g-coherent states coincide. 

3 Evolution equation of the g-nonentangled 
pure states 

In general, reduction of the pure state \ip > results in a mixed state p g and 
the unitary Schroedinger evolution of \ip > upon reduction becomes nonuni- 
tary, resulting in different forms (under different approximations) of master 
equations for p g (t). However, if the g-nonentangled pure state \ip > evolves in 
the subset of g-nonentangled states the reduced state always remains pure. 
In order for this to occur in general the Hamiltonian linear evolution of 
\ip{t) > is not enough and a nonlinear constrain has to be added to ensure 
the preservation of the g-purity P g (ip(t)). In order to formulate such con- 
strained evolution we shall use the classical Hamiltonian formulation of the 
quantum evolution. 

It is well known (please see [12] or [T5] and references therein) that the 
evolution of a quantum pure state in H N as given by the Schroedinger equa- 
tion can be equivalently described by a Hamiltonian dynamical system on 
R 2N with the evolution equations in the Hamiltonian form: 

x l = u^VjU, (4) 

where x i = q l = (c* + C;)/V2, % = 1, 2 . . . JV; x i = p l = ^l{c* - Ci)/y/2, i = 
N + 1, 2 . . . 2N is the vector of coordinates qi and momenta pi, and q are 
complex expansion coefficients of the pure state \ip > in some basis. The 
Hamilton's function "H(x) is given by the quantum expectation of the Hamil- 
tonian H in the state \ip >: H =< ip\H\ip >, and the inverse of the symplectic 
form uj % i is given by the imaginary part of the scalar product in H N . In the 
canonical coordinates Xj the symplectic form Uij assumes the standard form 

where and 1 are N dimensional zero and unit matrices. 

We shall use the classical geometric formulation of a quantum dynamical 
system in order to derive an equation for the quantum evolution constrained 
on a submanifold of R 2N that corresponds to pure coarse-grained states. 

Consider first the example of a pair of qubits. In this case the subspace 
of product states l^ 1 > ®\ip 2 > is characterized by the following condition: 



C1C4 = C2C3 on the expansion coefficients in the computational basis, which 
can be translated into a condition on the real canonical coordinates. This 
condition characterizes the submanifold of R 2N of points corresponding to the 
product states. The characterization of product, i.e. nonentangled, states is 
also given by the condition of maximality of P g (ip) where g is the algebra of 
local operators generated by o\ yz <g> 1, 1 <g> <y\ yz - In terms of the canonical 
coordinates (q,p) we obtained the purity function P g (q,p). In fact local 
purities < a]: 2 > 2 + < o^ 2 > 2 + < o^' 2 >2 as functions of the canonical 
coordinates are equal and are represented by half of the total purity function. 
It can be demonstrated that the function P g (p, q) is maximized precisely when 
the separability constraint C1C4 = c 2 c 3 is satisfied 

Hamiltonian equations with the algebraic constraint C\C^ = C2C3 have 
been studied for the first time in [13] and in [H]. The formalism of quan- 
tum constrained dynamics developed in [13] is based on Dirac's approach 
to classical constrained Hamiltonian systems and requires the constraint to 
be given explicitly in terms of an even number of independent real func- 
tions. In our present case there is in general only one real constraint: 
$(x) = P g {x) — max = 0, and the approach with symplectic constraints 
of [13], [H] can not be applied. However, the formalism of the so called 
metrical constraints, developed in [15], with an arbitrary number of real 
functions defining the constraint submanifold is applicable. For the benefit 
of the reader we shall rederive the constrained dynamical equations with only 
one real constraint which is of interest here. 

The purity constraint 

$(x) = P g (x) — max = (6) 

represents a single scalar condition that we want to impose on the evolution. 
In order to impose this condition the component of the Hamiltonian vector 
field x (4) normal to the constraint submanifold has to be removed resulting 
in 

x i = uj^WjH - Xg^Vj®, (7) 

where g^ is the unit metric on R 2N and A is a single Lagrange multiplier to 
be determined. Substitution of (7) in $(a;(t)) results in 

uj^Vi^VjH = A^'V^V^, (8) 

from which 

g^Vi^Vj® ' 

Substituting this A in (7) results in the constrained dynamical equations 

i^^H- ^Z'lZ 1 * ^!*- (10) 



We propose the reduction of the constrained equation (10) on the con- 
strained manifold to represents dynamical equation of the coarse-grained pure 
states. 

Observe that the numerator in (10) represent the Poison bracket {$, %} = 
$ and the denominator is ||V$|| 2 . Using the equalities 

T -A 9<L> T -A 9<L> m\ 

^ijQj — °ij g~ > ^ijPj — °ij g~ > K 11 ) 

where Lij are matrix elements of the operator L and in our case §(ip) = 
P(ip) — max the denominator can be further transformed as follows 



g" J2Vi<Li> 2 V j <L k >' 



2 V7 ^ T ^2 

l,k 

4J2< UL k >-<L l xL k >= 4X)(AL,) 2 = 4A(^). (12) 

l,k I 



Before presenting few examples we would like to make some comments 
concerning the constrained equation (10). 

1° In the open system picture of the distinguished system, to be discussed 
in the next section, and in the usual weak coupling approximation (WCA), 
with the distinguished observables identified with the Lindblad generators, 
the above equation is greatly simplified. Namely, in the WCA the Hamilto- 
nian and the system operators Li that couple with the environment operators 
satisfy 

[H,U] = X l L l (13) 

and $ satisfies 

* = PM=2£(AL,) 2 = 2AM, (14) 

so that the equation (10) is reduced to 



x i = (J'VjU - -zg^Vj®. (15) 



The open system interpretation of the constrained equation (10) and the 
equation (15) will be discussed in more details later in the next section. 

2° An equivalent constrained equations, of the same form as (10) and in 
the special case (15), are obtained if instead of the purity constraint P(i[>) = 
max the constraint A(^) = X)z(ALz) 2 = min is used. In particular, the 
special case equation (15), valid under the same conditions, with the use of 
(11) can be written in the form 

< ^- = -i[H,4>] + £(L 2 + < U > 2 -2<L l >L)\^>. (16) 



3° Number of variables and equations in (10), and in (15), can be reduced 
if, prior to imposing the constraints, the normalization of < ip\ip > and 
the global phase invariance of \ip > are explicitly used. The Hamiltonian 
Schroedinger equation (4) is then formulated on S 2N ~ 1 /S 1 instead of R 2N . 
The constrained equations have the same form as in (10) with the appropriate 
symplectic uo 1 ^ and metric g % i forms. An example is provided in the example 
b) below. 

4° The geometric Hamiltonian formulation and the constrained equations 
can be generalized to an infinite dimensional Hilbert space. 

Before we analyze an open system physical model of the coarse-grained 
dynamics let us present few examples of the g-constrained systems. 

Examples 

a) The first example is trivial in the sense that all pure states are g- 
coherent, and serves the purpose of illustrating the self-consistency of the 
approach. Consider a single qubit with the Hilbert space H 2 and an arbitrary 
Hamiltonian H . As the algebra of distinguished observables we take g = 
su(2). The g-purity is P g (ip) =< a 2 > + < a 2 > + < a 2 > and is maximal 
for any pure state. In this case all pure states are g-nonentangled and g- 
coherent. 

The constraint $ = P g {ip) — max = in the real canonical coordinates 
assumes the following form 

{p\+pl + ql + ql) 2 = 2. (17) 

The gradient of the constrains is given by 



v gi $ = 


= qM + pl + ql + ql) 


v g2 $ = 


= qiipl + pi + ql + q*) 


V P1 $ = 


= qi{p\ + pl + qi + ql) 


v P2 $ = 


( 2 , 2 , 2 , 2\ 

= P2{P! + P 2 + q! + q 2 ) 



and the Poisson bracket {$, %} is zero for arbitrary Hamiltonian % =< H >. 
Thus, the constraints are trivially satisfied, and the constraint dynamics is 
reduced to the linear Schroedinger part: x % = u^VjH. This example extends 
to the general case of H N with the distinguished algebra g = u(N). 

b) As the second example we consider the system of two qubits with 
the distinguished algebra of local observables g = su(2) ® su(2). In this 
case g-entanglement is the standard bipartite entanglement, g-nonentangled 
are the product states. Subsequent formulas are simplified if the condition 
< if}\ip >= 1 and the phase invariance are explicitly used. With this the 
system is reduced on the projective space S 7 /S 1 . Purity P(ip) =< (<r*) > 2 



+ < K) > 2 + < W) > 2 + < {<rl) > 2 + < fa 2 ) > 2 + < fa 2 ) > 2 in 
the computational basis and in canonical coordinates {qi, q2, <l3, Pi, P2, P3} of 
S 7 / S 1 is represented by 

P{q,p) = l + 4(2V2 Pl (p 2 q 3 -p 3 q 2 ) + 2q 2 3 + P 2 1 ( P 2 2 + q 2 2 )) 

+ 4{2V2 qi (p 2 p 3 + q 2 q 3 ) - q\{p\ + q\) - 2p\) (19) 

Pfa) = max is equivalent to C1C4 = c 2 c 3 where c\, c 2 ,c 3 ,c^ are coefficients 
of \ip > in the computational basis, and the equation of this constraint is 
equivalent to two real equations 

\/2]>i = P2<?i + Pi 92, v 7 ^ = q\q 2 - P\P 2 . (20) 

As for the Hamiltonian we consider two typical examples 

H, = a\ + a 2 z + nvlvl (21) 

H ns = a\ + al + fialal (22) 

(23) 

The reduced g-constrained dynamics in -fa, q 2 , pi,p 2 } coordinates is equiv- 
alently described by the equations (10) and the metrical constraints (19), or 
by the symplectic constrained equations with the constraints (20). The con- 
strained equations with constraints (20) turn out to be of a simpler form and 
are reproduced here. The details of the derivation have been presented in 

H3I- 

For the Hamiltonian H s the constrained equations read 



4:fip iqi q 2 + 2a;pi[2 + {p 2 f + fa) 2 ] 

2 + (P2) 2 + fa) 2 

4:fip 2 qiq 2 - 2cop 2 [2 + ( Pl ) 2 + fa) 2 ] 

2 + fa) 2 + fa) 2 

2/«7 2 [fa) 2 - (Pi) 2 - 2] + 2u; gi [2 + (p 2 ) 2 + (q 



<h 
Pi 

-/"/lH'fii — \l'21 — -! "I -^'/2|--i WW -r WW ! , 

P2 - , jz \2 , — \ 2 • y ZA > 



2) 



2 + (p 2 ) 2 + fa) 2 
2/igi[(g 2 ) 2 - (p 2 ) 2 -2] + 2a;g 2 [2 + ( Pl ) 2 + fa) 21 
2 + (pi) 2 + fa 



and for the Hamiltonian H ns . 

2^i[(p 2 ) 2 + fa) 2 - 2)] - 2w Pl [2 + (p 2 ) 2 + fa) 



qi 
q.2 



2+(p 2 ) 2 + fa) 2 
2m[(Pi) 2 + (qi) 2 - 2] - 2 W p 2 [(2 + ( Pl ) 2 + ( qi ) 2 } 



2+( Pl ) 2 + fa) 2 
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-2 Ml [(q 2 ) 2 + (p 2 f - 2] + 2a; gl [2 + (p 2 ) 2 + (g 2 ) 2 ] 
Pl ~ 2 + (p 2 ) 2 + (g 2 ) 2 

-2/xg 2 [fa) 2 + ( Pl ) 2 - 2] + 2coq 2 [2 + ( Pl ) 2 + fa) 2 ] 
P2 " 2 + ( Pl ) 2 + te) 2 • (25) 

There are also the equations expressing q 3 and p 3 in terms of qi,q2,Pi,P2, 
but the solutions of these are already given by the constraints. 

The dynamics generated by (24) and (25) is illustrated in figures 1 and 
2. In fig. 1 we illustrate the time series qi(t) for single typical orbit of 
(24) (fig. la) and of (25) (fig. 1 b). In figures 2a,b,c,d the Poincare sections 
q 2 = 0,p 2 > for H ns (25) are shown. It should be observed that g- 
constrained dynamics of the symmetric Hamiltonian H s is regular, while that 
of the Hamiltonian H ns with no such symmetry displays typical properties 
of the Hamiltonian chaos. Thus, although the linear Schroedinger equation 
always generates an integrable Hamiltonian system, the coarse-grained quan- 
tum system evolving according to the constrained equations can display all 
complexities of typical chaotic dynamics. 

c) In this example we again consider a system with g = su(2) distin- 
guished algebra but with the spin s = 1 i.e. with H 3 Hilbert space. As for 
the Hamiltonian we take a nonlinear expression of su(2) generators 

H = J z - 2 J x + fxJ 2 (26) 

When /i ^ the Schroedinger evolution with the Hamiltonian (26) does 
not preserve the sw(2)-coherent states. The set of sw(2)-coherent states is 
preserved when \x = 0. 

The g-constraint $(V>) = P su ( 2 )(^) - 1 =< J 2 > + < J 2 > + < J 2 > -1 
in the eigenbases of J z and in the real canonical coordinates of R 6 assumes 
the form 

^Psu{2) (q, p) = -4 + p\ + p\ - 2p\q\ + q\ + 8p 2 p 3 q ± q 2 + 2p\q\ 

+ 2q\q\ + 2p 2 2 {pl + ( Ql - q 3 ) 2 ) + 4g ig 2 g 3 + 2{pj - q\ + q\)ql 

+ <?s + 4 Pi (P2P3 - Psql + 2p 2 q 2 q 3 ) 

+ 2pi(pi-p3 2 + q 2 1 +q 2 2 -ql) (27) 

The Poison bracket of the constraint and the Hamiltonian, that is needed 
for the constraint equations (10), reads 

(J'V&VjH = 2fi[{p 3qi + piq 3 )(ql - p\) + 2p 2 q 2 {p lP3 - q lQ3 )]. (28) 

We see that, when \x = 0, the Poison bracket (28) is zero and the g- 
constrained equations reduce to the Schroedinger equation. The squared 

9 



norm of the $ gradient is given by somewhat complicated function of the 
canonical coordinates (q,p) and will not be reproduced here. We illustrate 
the form of simplified constrained equations (15) by the formula for p\ 

P! = -V2[(l + n) qi -q 2 ] 

+ Piipl -Pi + ql + <ll - 5s) - 2^292^3 + Psql ~ p\pz ~ Pi, (29) 

where the first line is the Hamiltonian term and the second line is from the 
gradient of the constraint. 

We shall come back to this example in the next section. 

4 Open system model 

The coarse-grained system specified by the distinguished variables can be 
considered as an open system with the larger closed system characterized by 
the full algebra u(N). In the case when the Hilbert space can be split into the 
tensor product with one component corresponding to the distinguished re- 
duced system the standard open system model of decoherence applies. This 
theory singles out a distinguished set of states, the pointer or the robust 
states, and characterizes them as pure states of the reduced open system that 
remain pure under evolution, or as states in which the reduced open system is 
not and does not get entangled with the environment during the full system 
evolution. Reduced states of the general coarse-grained system, discussed 
in sections 2 and 3, satisfy the same properties as the robust states of an 
open system under decoherence if the interaction of the open system and the 
environment is mediated by all of the distinguished observables. It has been 
demonstrated that the robust states in this case coincide with the g-coherent 
states [TU] in the weak coupling limit. In this picture the coarse-graining 
physically occurs due to decoherence of the distinguished system induced by 
specific interaction with the environment which represents generalized simul- 
taneous measurement of all distinguished observables. The pointer states are 
identified with reduction of the g-nonentangled or g-coherent states. 

We would like to identify the distinguished observables with observables 
that are simultaneously measured on the open system. In the weak coupling 
limit (WCL) the Born-Markov and rotating wave approximations result in 
the Lindblad master equation of the open system dynamics [16J 



pit) = -i[H, p\ + W ([Up, L\\ + [L hP L\) , (30) 



2 , 



where H is the open system Hamiltonian and L\ are the so called Lindblad 
operators. L\ are the open system operators that are coupled with that what 
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is considered as environment. If the eq (30) corresponds to the measurement 
of certain observables than Li are the Hermitian operators that represent 
the measured observables. In our model of the coarse-graining we shall sup- 
pose that the distinguished algebra is precisely the algebra formed by the 
Hermitian Lindblad operators in (30). 

As pointed out, the pointer or robust states in the open system model of 
decoherence are the pure states of the open system that remain pure in the 
course of evolution. It has been suggested j9] that an approximate evolution 
equation of the pure robust states can be obtained by minimizing the Hilbert- 
Schmidt (HS) distance from p(t) to the set of pure states. In the case p(t) is 
given by the Lindblad eq. with Hermitian Lindblad operators the equation 
of the HS closest pure state is [H] 

^Ll = -i[H,i>} + £(L?+ < U > 2 -2<L t > L«)|V > • (31) 



This is precisely our constrained equation (15) when g-entanglement measure 
P(ip) = J2i < Li > 2 is replaced by the equivalent measure A(ip) = J]/(AL;) 2 
in the WCL with the distingushed observables being the Lindblad generators. 
The equation (31) (or (15)) represent the deterministic part of the stochas- 
tic Schroedinger equation derived in the quantum state diffusion theory [T7] 
for arbitrary random pure state, which we reproduce here because it will be 
used for numerical computations in the next example. The Ito form of the 
QSD equation corresponding to (30) reads 



\dip > = -iH\ip > dt 



+ 



£2 < l\ > u - l\u- < l\ >< u > 



L / 



\if)(t) > dt 



+ J2( L i-< L i>)\' l P( t )> dW i ( 32 ) 

1 

where dWi are independent increments (indexed by I) of complex Wiener 
c-number processes Wi(t) satisfying 



[dWt] -- 


= E[dWidW v ] 


= 0, 




dWidWv 


= $i,i>dt, 




I 


= 1,2. ..m 



(33) 

where E[] denotes the expectation with respect to the probability distribu- 
tion given by the (m-dimensional) process W, and W\ is the complex conju- 
gate of W\. 
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The random vector \ip(t) > which satisfies (32) is related to the density 
matrix p(t) which satisfies the Lindblad equation (30) by averaging over the 
realizations of the process (32) 

P (t) = E[\m >< m\]- (34) 

Let us stress that the HS approximate robust state eq. (31) assumes 
validity of WCL and coincides with the constrained eq. (15) simplified from 
(10) under this assumption. On the other hand the general constrained 
evolution given by (10) is valid, in the sense that it preserves P(if>) and 
A(^), with no assumption about special evolution of P(ip) which is obtained 
under the WCL. 

An example 

In the case of an open system that satisfies the conditions for the weak 
coupling approximation (13) the dynamics of the system is described well 
by the Lindblad equation and the pointer states are exactly the g-coherent 
states [10] . Using particular examples, it has been demonstrated [9] that the 
equation (31), which coincides with the simplified form of the constrained 
equation (15), describes well the evolution of the pointer i.e. g-coherent 
states. We shall analyze here an example that does not satisfy the condition 
(13) of the WCA. 

Let us consider, as an example, the two mode Bose-Hubbard model (see 
for example [IS]), given by the following Hamiltonian with h = 1 

H = e\a\ai + e 2 a 2 a 2 + a(a[a 2 + a 2 ai) + p(a\ a 2 + a 2 a 2 2 \ (35) 

where ai,a\,i = 1,2 are bosonic annihilation and creation operators of the 
two modes. The dynamics preserves total particle number N = a\a\ + a 2 a 2 . 
Introducing operators 

Qj = (<4 + a i)/^i Pi = ^ a ) ~ a j)/^i J = !> 2 > ( 36 ) 

or the operators 

J x = -{a\a 2 + a\ai) 
J y = -(oi0 2 -a|oi) 

J z = -(a\a 2 -a|ai)), (37) 

the Hamiltonian assumes the following forms respectively in coordinates (36) 

H = t 1 {pl + ql)/2 + p{p\ + ql) 2 /A + 
+ e 2 {p 2 2 + q 2 2 )/2 + p{p 2 2 + q 2 2 ) 2 /4 + 
+ a(pip 2 + q x q 2 ) (38) 

12 



and in terms of (37) 

H=-2aJ x + 2{e 2 -e 1 )J z + pJ 2 . (39) 

In what follows we shall always set a, — 1, e 2 — e\ — 1. 

The preserved total number of particles is related to J 2 by J 2 = N/2(N/2+ 
1). Thus, the effective Hilbert space of the system carries an irreducible repre- 
sentation of SU(2), which is the dynamical group of the model. This suggest 
that the SU{2) coherent states have a special status in the model (35). This 
however is not true, because the nonlinear term pj 2 makes the set of SU(2) 
coherent states noninvariant. 

We would like to analyze system (35) interacting with an environment 
via operators (37) or (36). The Hamiltonian (35) and operators (36) or (37) 
used as the Lindblad operators do not quite satisfy the condition (13) for 
the WCA. Nevertheless, we shall suppose that the open system evolution is 
described by the Lindblad equation with Lindblad operators given either by 
(36) or by (37). Notice that the result A g (ip) —$■ min obtained in [TU] does 
not apply necessarily since the system does not satisfy the WCA condition. 
We shall demonstrate that the asymptotic states of the Lindblad eq. of an 
open BH system interacting with an environment via the Lindblad operators 
L[ satisfy the constraints condition A g (ip) = min almost exactly 

Let us first consider the open system evolution in terms of random pure 
states \^{t) > and the QSD equation (32). We first choose Li = J X ,L 2 = 
J y ,L 3 = J z and compute A su (2)(i()(t)) and P su (2){'4 ) (t)) from an initial state 
equal to the number state given by (ai) 2 (c4) 2 |0, >. The results are shown 
in figure 3. The state quickly converges to those with a minimal A su ^(ip(t)), 
i.e. to the sw(2)-coherent states. On the other hand A# 4 = A 2 pi + A 2 p 2 + 
A 2 qi + A 2 q 2 remains constant and large. sw(2)-purity is less than maximal 
at the beginning but quickly converges to the maximal value. Although the 
state \ip(t) > is always a pure state of the Hilbert space, its su(2)-purity is 
maximal only when \ip > is an sw(2)-coherent i.e. an su(2)-nonentangled 
state. Analogously, assuming the Lindblad operators to be L\ = qi,L 2 — 
q 2 ,L 3 = pi,L 4 = p 2 implies an evolution such that A su ( 2 ) is far away from 
its minimum, but A# 4 converges to values close to the minimal and remains 
such for almost all times (please see fig. 4). 

The equivalent conclusions are obtained when the evolution is described 
in terms of p(t) = E(\if)(t) >< ip(t)\), i.e. by the Lindblad equation. This 
is illustrated in figure 5 and 6, with A g for g = su{2);g = H^ and L\ = 
J xi L 2 = J y , L 3 = J z with the number initial state (fig. 5), and L\ = qi, L 2 — 
q 2 ,L — 3 = pi,L 4 = p 2 with an su{2) coherent initial state in figure 6. 
Only two hundred QSD sample paths are use to compute pit) and then the 
corresponding A g (p) 
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Furthermore, consider evolution from an su(2) coherent initial state with 
J = 1 and with the Lindblads being J x , y , z - The Lindblad eq. assumes WCL, 
the asymptotic states satisfy A su ^) ~ min, and the simplified constrained 
equation (15) applies. Indeed, the Lindblad evolution is well approximated 
by the simple form of the constrained equation (15), as is illustrated in figure 
7a,b. 

The usual picture of decoherence applies: An arbitrary initial state evolves 
very quickly into a mixture of g-coherent states, and then each of these 
evolves in a way that is well approximated by the nonlinear eq. (10) or in 
the WCL by (15). 

Let us stress that the coarse-graining by distinguished observables, dis- 
cussed here, is specially appropriate in a description of macroscopic features 
of a quantum system, with the distinguished observables identified with the 
macroscopic quantities. In this case the Hilbert space of the quantum sys- 
tem does not have the bipartite tensor product structure, with one party 
being characterized by the macroscopic observables and the other party be- 
ing the environment. The usual models of decoherence [Tj5] with the initial 
separation \ip >= \ip s > ®\ip env > do not apply. However, the picture of 
coarse-graining by distinguished observables with the corresponding nonlin- 
ear evolution equations can be applied. 

5 Summary 

We have analyzed the coarse-graining introduced by a chosen set of distin- 
guished observables. The algebra of distinguished observables defines the 
corresponding generalized nonentangled states which coincide by definition 
with the generalized coherent states. The states obtained by reduction on 
the distinguished observables of the g-nonentangled states are pure. We have 
propose to consider the coarse-grained evolution as constrained Schroedinger 
dynamics, where the constraint guaranties that the state is always pure g- 
nonentangled. In order to formulate the constrained evolution equations we 
used Hamiltonian formulation with the metrical form of the constrained dy- 
namics as developed in [TB] . 

Further on we discussed an open system model of the coarse-graining and 
of the reduced constrained equation. In the weak coupling limit the open 
system dynamics is given by the Lindblad master equation. In this limit, 
and if the Lindblad operators are taken to represent the distinguished ob- 
servables then the constrained equations for the g-coherent states developed 
here coincide with previously suggested [9] evolution equation for the pointer 
states of the open system. 
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Our simplified constrained evolution equation (15) for the g-nonentangled 
states coincides with the deterministic part of the Ito stochastic Schroedinger 
equation developed in the quantum state diffusion theory (QSD) of open sys- 
tem dynamics. The stochastic Schroedinger equation describes dynamics of 
any random pure state. Our constrained equations describe dynamics of de- 
terministic pure states which are in the subset of all pure states that remain 
pure during the evolution. From a formal point of view, it would be interest- 
ing to derive the QSD stochastic equations using the formalism of constraints, 
where the constraint would be given by random variables representing the 
obtained results of measurements with Gaussian distribution. 

The formalism of coarse-graining as the constrained evolution can be used 
to study coarse-grained macroscopic observables of a quantum system and 
derive their classical behavior. 
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FIGURE CAPTION 

Figure 1. Illustrated are segments of time series qi(t) for the system (24) 
in a) and (25) in b). 

Figure 2. Poincare sections q 2 = 0,p 2 > and % = 1.5 for the system 
(25). The parameters are (a)/i = 1.1, (b) /i — 1.3, (c) fi = 1.5 and (d) /i — 1.7 

Figure 3. Illustrates the invariant fluctuation A s (^) (2) in the cases 
g = sw(2) (full line) and (7 = H4 (dotted line) for the QSD evolution with 
the Hamiltonian (35) and Lindblads L\ = J x , L 2 = J y , L% = J z . The initial 
state is the number state |2,2 >= (g4) 2 (g4) 2 |0, >. The parameters are 
H = 0.l,a = l,e = 0,7 = 0.9 

Figure 4. Illustrates the invariant fluctuation A s (■?/>) (2) in the cases 
g = su(2) (full line) and g = H^ (dotted line) for the QSD evolution with 
the Hamiltonian (35) and Lindblads L\ = q l7 L 2 = q 2 , L 3 = p l7 L 4 = p 2 . The 
initial states is an su(2) coherent state. The parameters are \x = 0.1, a = 
l, e = 0,7 = 0.9 

Figure 5. Illustrates the invariant fluctuation A 9 (p) (2) in the cases g = 
su(2) (full line) and g — H 4 (dotted line) for the evolution by the Lindblad 
eq. with the Hamiltonian (35) and Lindblads L\ = J x , L 2 = J y , L% = J z . The 
initial state is the number state |2,2 >= (a{) 2 (o4) 2 |0, >. The parameters 
are /i — 0.1, a — 1, e = 0, 7 = 0.9. 

Figure 6. Illustrates the invariant fluctuation A 3 (p) (2) in the cases g = 
su(2) (full line) and g = iJ 4 (dotted line) for the evolution by the Lindblad eq. 
with the Hamiltonian (35) and Lindblads L\ = qi,L 2 = 52,-^3 = Pi, L4 = p 2 . 
The initial states is an su(2) coherent state. The parameters are \x = 0.1, a = 
1,6 = 0,7 = 0.9 

Figure 7. Evolution of Tr[pa z ] according to the Lindblad eq. (30) (full 
line) and of < a z > according to the simplified constrained eq. (15) (dotted 
line) with the Hamiltonian (35) and Lindblads L\ = J X ,L 2 = J y ,L 3 = J z . 
The initial states is su(2) coherent state \j,j z >= |1, — 1 >. The parameters 
are // = 0.1, a — 1, 7 = 0.2 and a) e = and b) e = 1. 
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